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Abstract. Due to the increase in the poloidal wave momentum of the lower hybrid 
wave as it propagates into the tokamak, the contribution to the electron toroidal 
angular momentum of the parallel momentum resonant transfer is larger than the 
initial toroidal angular wave momentum. A significant contribution in the opposite 
direction is given by the perpendicular momentum transfer that makes the total angular 
momentum transfer equal to the initial toroidal angular momentum, and it results in 
the outward radial pinch of electrons 
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1. Introduction 

The momentum of radio-frequency (RF) waves has been studied since the early days of 
development of electromagnetics pQ. Recently, the experimental observation of plasma 
flows generated by RF waves has renewed the interest in momentum deposition by RF 
waves. For example, a significant ion toroidal rotation (~50km/s) has been measured 
by X-Ray spectroscopy for impurities in Alcator C-Mod during lower hybrid (LH) wave 
power injection [2], and the relation between the computed toroidal angular momentum 
input from LH waves and the measured initial change of ion toroidal rotation has been 
investigated [3]. In tokamaks, the lower hybrid wave is used to drive plasma current 
with asymmetric antenna spectra along the direction parallel to the static magnetic field 
[I] . Due to electron Landau damping, the wave power is transferred to the non-thermal 
fast electrons (vn ~ 3v te — 10t> te , where vn is the parallel electron velocity and v te is the 
electron thermal velocity) [U EJ |6] . The toroidal phase velocity of the wave is chosen 
to increase current drive efficiency and ensure the accessibility of the wave to the core 
of the tokamak. The electron Landau damping can be described kinetically as a quasi- 
linear velocity diffusion coefficient if the strength and the spectrum of the electric field 
satisfy some conditions given in EJ |9]. 

In this paper, we investigate how the wave momentum changes as the wave propa- 
gates from the launcher to the core of the tokamak where it is damped and transferred 
to the plasma. We also investigate how the transferred momentum affects the radial 
motion of the particles. Studying the wave momentum transfer is essential to explain 
the temporal behavior of the ion toroidal rotation initiated by the LH wave injection. 
The turbulent momentum transport dominates the temporal evolution on the transport 
time scale [3j [10] , but before studying this long time scale behavior we need to under- 
stand the momentum deposition. 

Figure 1 shows the typical behavior of a LH wave in an inhomogeneous tokamak. As the 
wave propagates from the low field side launcher, it develops a very high poloidal wave 
vector (more than 10 times larger than toroidal wave vector) due to the plasma disper- 
sion relation [TT]. The large poloidal wave vector contributes to the parallel wave num- 
ber k\\ as much as the toroidal wave vector does, even overcoming the small ratio of the 
poloidal magnetic field over the toroidal mag netic field, B0/B4 ~ 0.1 (k\\ ~ ^k^ + ^-kg 
in a circular tokamak, where B^, Bg, B are toroidal, poloidal, and total magnetic field, 
and k^, kg are toroidal and poloidal wave numbers, respectively). That results in the 
parallel refractive index nil = h\c/oj ~ —3 of the damped wave, significantly larger than 
the toroidal index = k^c/u ~ —1.9 at the launcher, as shown in Fig. 1 (the negative 
sign means that the wave propagates in the counter-current direction of the tokamak). 
The frequency u is the wave frequency, and c is the speed of the light. The electron 
Landau damping of the wave becomes stronger where the phase velocity of the wave 
becomes lower (in other words, where the refractive index becomes higher), since a lower 
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phase velocity resonates with more electrons. 

As shown in Fig.l, until the wave reaches the region where the parallel phase velocity of 
the wave is sufficiently reduced by the poloidal coupling (e.g. n\\ ~ —3) to interact with 
less energetic electrons, the resonant interaction is negligible. Nevertheless, the poloidal 
momentum of the wave changes due to the inhomogeneity of the magnetic field and the 
plasma density and temperature. There is a significant poloidal wave momentum gain. 
The wave gains poloidal momentum slowly in the non-resonant region, and then trans- 
fers it in a short distance where it resonates. However, the toroidal angular momentum 
of the wave does not change due to the toroidal symmetry, and the original amount is 
fully transferred to the plasma in the resonance region (see the constancy of the green 
line in Fig.l). The non-resonant interaction can be studied as a combination of the 
Reynolds stress and the Lorentz force in both fluid models [121 13] and kinetic models 
[HI [15j [161 [17]. The equivalent non- resonant ponderomotive force has been evaluated 
in different ways, and it has no effect on the toroidal flow (T7J HE] ■ 

When the wave energy is transferred to the the plasma due to a resonance, the cor- 
responding wave momentum (k/u times the energy density, where k is the wave vector) 
is also transferred to the plasma [1]. This relation has been verified by evaluating 
the Lorentz force in fluid models [T2l IT9] and kinetic models [20| ETJ [22] . However, 
the momentum transfer by resonance has been calculated incorrectly for the LH wave 
[31 [231 El] resulting in an incorrect radial electric field. These calculations have ignored 
an important contribution to the Kennel- Engelmann quasilinear diffusion coefficient. 
The Kennel-Engelmann quasilinear diffusion coefficient [5] describes the resonant in- 
teraction of the plasma with the wave. The gyroaverage of this quasilinear operator 
is used to model the diffusion of the distribution function in velocity space. However, 
since some components of the momentum, such as the toroidal direction, depend on the 
gyro-phase, the diffusion in gyro-phase must be taken into account for momentum trans- 
fer calculations, and the gyroaveraged quasilinear operator is not sufficient to explain 
the total momentum transfer. We reexamine the amount of momentum transfer from 
LH wave to the plasma by resonant interaction in this paper. The new contribution 
to momentum transfer that we find is important because the poloidal wave number is 
large in the resonance region, giving k\\ > k^ as we have discussed above. Using the 
gyroaveraged quasilinear operator only transfers the parallel wave momentum, leading 
to an incorrect evaluation of the toroidal angular momentum transferred by the wave. 

Once the momentum of the LH wave is transferred to the electrons, most of it is trans- 
mitted to the ions by electron-ion collisions, but the rest is balanced by an electron 
radial pinch. The radial non-ambipolar electron pinch has been proposed as an expla- 
nation for the ion rotation induced by LH waves [21 [251 126]. It has been argued that 
the counter-current direction momentum transfer from the LH wave to the trapped 
electrons induces a radially inward pinch, and it results in an additional inward radial 
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toroidal angle[7i rad] 



Figure 1. Evolution of the toroidal angular momentum of the LH wave in terms of the 
propagation toroidal angle from the launcher. The solid line is the toroidal projection 
of the parallel refractive index multiplied by the major radius (n\\(B$/ B)R), which 
is an important parameter for Landau damping. The color of the solid line is the 
normalized Poynting flux of a ray. The power of the LH wave is absorbed by electron 
Landau damping beyond toroidal angle 3.57T where its color changes from red to blue. 
In this power absorption region, Tin (B^/ B)R is around -2.3, much higher than original 
toroidal refractive index multiplied by the major radius (n^R ~ —1.4). This graph 
corresponds to one of the LH wave rays for Alcator C-Mod with B^ — 5.3T, plasma 
current I p = 700fcA, major radius Rq = 0.67m, minor radius a = 0.22m, electron 
temperature T e (0) = 3.5KeV, electron density n e (0) = 1.2 x 10 20 , initial nil = —1.9 
and P a bs = 0.8MW. These profiles are calculated with Genray-CQL3D [28 . 

electric field to ensure ambipolarity. This excess radial electric field will then act on the 
ions, leading to ion rotation [25]. In this paper, we propose another (actually stronger) 
mechanism that gives an outward electron pinch. 

The rest of this paper is organized as follows. In Sec. |2} we revisit the quasilinear 
diffusion operator taking into account the flux in the gyro-phase angle direction to eval- 
uate the total amount of momentum transfer by resonant particles. In Sec. [3j we discuss 
briefly the wave momentum gain or loss by non-resonant effects. The inhomogeneity of 
the tokamak system generates poloidal wave momentum. In Sec. HJ we discuss the radial 
non-ambipolar electron pinch due to the resonant momentum transfer. The pinches of 
the passing electrons and the trapped electrons due to the LH wave parallel momentum 
are explained. More importantly, we present a new mechanism due to the perpendicular 
momentum transfer that gives an outward electron pinch. Finally, the conclusions of 
this paper are given in Sec. El 
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2. Momentum transfer by resonance 

In this section, we reconsider the quasilinear diffusion in velocity space including the 
gyro-phase. We show that we must retain the contribution to the momentum transfer 
from the quasilinear diffusion in the gyro-phase direction. In this section, we only focus 
on the momentum transfer by resonant particles and its relation to the power absorp- 
tion. We derive a quasilinear diffusion operator, and we use it to evaluate the momentum 
gained by the particles by applying the resonance condition. This proof can be applied to 
any type of resonance (cyclotron or Landau damping) and any direction of momentum. 
For convenience and without loss of generality, we discuss toroidal angular momentum 
which has both perpendicular and parallel components. Let z be the direction parallel 
to the static magnetic field, and x, y the orthogonal coordinates (see Figure 2). Then, 
using the gyro-phase angle a, the velocity is v = v± cos ait+v± sin ay+vnz = Vxp+Vnz, 
and its toroidal component is = v± coso;(x ■ <f>) + v± sina(y ■ <j>) + i>ii(z ■ </>), where 
p = v±/v± and </> is the unit vector in the toroidal direction. The wavenumber vector 
is defined as k = k± cos /3x+ k±_ sin /3y + k»z — k± cos (a — (3)p — k± sin (a — (3)a + fciiz, 
and the electric field is E = E x ~k + E y y + E\\z = E±p + E a a + E\\z, where a = 
z x p is the unit vector perpendicular to both vi and the magnetic field. Here, 
E± = E x cosa + EySma = (E + + EJ) cos (a — 0) — i(E + — EJ) sin (a — /3), and 
E a = —E x sin a + E y cos a = —i(E + — EJ) cos (a — /3) — (E + + EJ) sin (a — where 



where the triangular bracket (...) in (TjQ) indicates the average over a number of wave 
periods in time and space. Here, m and Ze are the mass and the charge of the species of 
interest, respectively, e is the charge of the proton, and I is the unit tensor. We have used 
the Fourier analyzed perturbed fluctuating electric field, E = Ek exp(ik-r— iu-kt), the 
fluctuating magnetic field B = ^ k Bkexp(zk-r — iu^t), and the fluctuating distribution 
function, / = ^ k / k exp(zk • r — ko-^t). The functions E k = E(a> k , k), B k = B(w k ,k), 
and /k = /(a; k , k) satisfy the relation /_ k = /(a;_ k , — k) = /*(a; k ,k) where * denotes 
complex conjugate and = — w* k - Faraday's law has been used in going from (JTJ) to 



\{E x ±iE y )e^. 



The quasilinear diffusion operator is obtained from 




(1) 




6 



to write B k = (c/a>)k x E k f. The quasilinear operator can be written as 

1 d 



Q(f) - - 

m 



v± ov± 

The flux in the perpendicular direction is 

k\\V\\ 
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1 dT a dY\ 
H - 

v±_ da dv\\ 

k±v 



K,\\— — cos (a - /3))/ k , 

LU J 



(3) 



(4) 



the flux in the gyro-phase direction is 



T a = -^{£* )Q (l-^cos(a-/3)- 



h v \\\ 



UJ 



El 



k±v± . 



and the flux in the parallel direction is 

k ± v± 



-EK.il 1 



sin (a - 0) - £ k „ sin (a - 0) }/ k 



cos (a - f3) ) + 



(5) 



(6) 



Here, the perturbed fluctuating distribution function consistent with a single mode wave 
[31] is 



Ze f 

/k = exp(— zk ■ r + £a>t) / dt' expfik ■ r — iut') 

m J-oo 

V«'/o, 



x E k - 



I I 



v'k 



-oo 

v'k 



(7) 

where (£', r', v') is a point of phase space along the zero-order trajectory. Its end point 
corresponds to (t, r, v). The background distribution, / = fo(t,r,v±,v\\), is gyro-phase 
independent because of the fast gyro-motion. As a result, 

Ze f°° r 

/k = / drexp(2 7 ) cos (77 + fir)((£ kj+ + E K _)U - ^V) 

m Jo ^ 



1 sin (77 + Qr)(E K+ - E^)U + ^k,||^}- 



(8) 



Here, r = t—t', and 7 = (a;— \\)t— A(sin (77 + fir)— sin (77)), where A = -^fS 77 = a—/3, 
Q = ZeBo/mc is the gyrofrequency and Bo is the magnitude of the background mag- 
netic field. Also, U = |k + ^ (Vl£& - unl^Y and V = ^ (v ± $& - uiil^V We 

ofx w \ oily N ovx / ' cj \ - 1 - c«;|| II iroi 1 

follow Stix' notation in 1311. 



f In typical tokamak geometry, the toroidal and poloidal spectra are discrete due to periodicity, but 
the radial spectrum is continuous. Also, the parallel spectrum is continuous for the flux surfaces 
with non-rational safety factor by the coupling of toroidal and poloidal components [9j. Even though 
using integrals in Fourier space would be more appropriate, we use the notation ^ k for simplicity. 
The summation in the discrete toroidal and poloidal spectrum space is also closer to the numerical 
evaluation in a code [29j [30] . 
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For the energy transfer, the contribution of the flux in the gyro-phase direction vanishes 
due to the integral over a, 

,2 



abs 



dv\\ / dv±2nv± 

-oo JO 

OO /"OO 

dv\\ I dv±2nv± 

■oo Jo JO 



mv 



W) 



277 



da Zev 
2^~ 



2 r 



1 d , < dT 

(v±r ± ) + 



v± dvj 



dv\\ 



(9) 



with (...} a = ^ Jq* da{...) the gyroaverage. For this reason, the typical Kennel- 
Engelmann quasilinear diffusion operator |5] is gyroaveraged and does not retain the 
flux in the gyro-phase direction. For completeness, we have evaluated the energy trans- 



fer, P a bs, in Appendix A 



The gyroaveraged quasilinear operator is not sufficient to calculate the toroidal momen- 
tum transfer, which has gyro-phase dependent components. The total toroidal angular 
momentum deposited by the wave is 



/oo /*oo 
dv\\ / dv±.v ± (mRvfQ(f)) a = Pj + AP^ + AP^ 
-oo Jo 



where 



/■oo 

ZeR I dv\\ 

-oo 



2tt 



(10) 



(11) 



dv±v± / da(z ■ <f>)T\\ 
'o Jo 

is the component of momentum transfer that one obtains when using the gyroaveraged 
quasilinear operator, whereas 

c2tt 

APJ: = ZeR I dv\\ / dv±v± / da{— cosa(x • 0) — sina(y • 0))r_i_ , 

/*oo /*2"7r 

/ dv±v± / da(sma(St ■ 0) — cosa(y ■ 4>))T a 
'o Jo 

are the contributions that appear when the complete dependence on the gyro-phase 
is retained. Using the perturbed distribution function and the expansion in Bessel 



AP? 



-oo 

"OO 

ZeR I dvu 

-oo 



(12) 
(13) 



functions described in Appendix A , the toroidal momentum transfer term in the parallel 
direction, Pi, becomes 



nZ 2 e 2 R 



rn 



k 



/oo /*oo 
dv\\ / dv±2nv± 5(co — k\\v\\ — nf2)(z ■ 0) 
oo Jo 



XKn\ 2 Hfo) = {^Pabs >k R(Z ■ 0)} 



nil 



/ . P a bs,k R ( Z * 0) f ' 



(14) 



where Xk,n = E^\\J n -^ + Pk,+^n-i + P kr _J n+1 is the effective electric field, and 
J n (A) are the Bessel functions of the first kind with integer order n. The operator 
L(f ) = ( 1 - to) + dfsi i s introduced in [3TJ (see |AppendixT?1 for the 



detailed derivation). The piece of the momentum transfer Pj is directly related to the 
quasilinear diffusion operator used to calculate the power absorption (compare equation 
fl!4p with (IA.8p ). The direction of diffusion is determined by the characteristics of the 
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operator L(f ) in (IA.7I) (i.e. the tangents to the contours v\ + yu\\ — j^-J =constant), 
and the magnitude of the diffusion is determined by the projection of the distribution 
function gradient onto the these characteristics [51 [31] (see Fig 3). In particular, for 
electron Landau damping of the LH wave (i.e. u = k\\V\\), the piece of the toroidal 
momentum transfer p\ can be simplified to the following equation within a small error 

_ „2 D roo pod f\r 

pUld- E/ dv W dv ± 2™ L 8{u-k m )\E^^fi{\){l-fa 

m e k J-oo JO OV \\ 

= j°° dv\\ d Vl _2nv ± m e Rv^ (^(jjf^) a ^j , (15) 

where m e is the electron mass and the gyroaveraged quasilinear diffusion coefficient for 
electron Landau damping is (-D^ Lr> ) Q = ^-^2^5(0} ~ kll^lOI-^MP^oM- For Landau 
damping, the quasilinear diffusion happens only in the parallel direction (see Fig. 3). 
Equation ( fl5l) exemplifies the problems that appear if the gyroaveraged quasilinear 
diffusion operator is employed to evaluate momentum transfer. Using the typical 
gyroaveraged quasilinear diffusion coefficient, we can only evaluate the momentum 
transfer pjj instead of the full momentum transfer P^. For the rest of the toroidal 
momentum transfer, we need the quasilinear diffusion operator before the gyro-phase 
averaging, 

^ da 



/OO /»00 /» 

dv\\ / dv ± 2nv ± 
-oo JO Jo 



2n 
k\\V 



p* k ± v \\ 



UJ 



x (cos /3(x • 0) + sin (3(y ■ 0)) + (l - -^j (E^x • + E* k J ■ 0) 

+ % (^) (E* >+ - E^_) (sin a(x • 0) - cos a(y ■ 0))] / k . (16) 



Using Appendix B, we can simplify this equation to 



oo ^oo 



ap; + ap; = -^-^ E 1 dv \\ I dv ± 2 ™± 



oo JQ 



5(u - k {l v l{ - nfi)(cos/3(x • 0) + sin/3(y • 0)) -\xKn\ 2 L(fo) 

n 

= E (^T^ Pa6s ' kjR ) = E {^^ s , k P(cos/3(x • 0) + sin/3(y • 0))} . (17) 

In conclusion, for any resonance (e.g. cyclotron, Landau damping), the total toroidal 
angular momentum transfer according to f ll4p and (j!7p is 

p+=p\+ + ^ = E (^ p ^^) ' ( 18 ) 

as expected [TJ [121 US 1201 [211 122] . For electron Landau damping, this result implies that 
we need to consider the Lorentz force V\\z x Bj_ for momentum transfer, where Bj_ is 
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Figure 2. Sketch of a wave vector k = (m/r)6 + (n/R)<j> = kj_ + k|| in a parallel 
(z) -perpendicular (x, y) coordinate system and in a toroidal (^)-poloidal (9) coordinate 
system. For the LH wave, the component fcii has a bigger toroidal projection than the 
initial toroidal component n/R at the launcher due to the poloidal coupling. This 
difference is balanced by the toroidal projection of k± which has an opposite toroidal 
projection to the original n/i?. 



induced by Faraday's law, as shown in Fig. 4. This mechanism supplements the parallel 
direction transfer by the parallel electric field En to give the whole momentum transfer. 
In other words, the Lorentz force in the perpendicular direction described by E£ u 
in (flBjl results in the perpendicular acceleration that we need to satisfy equation (1TB]) . 
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<Sv kJ > a =0 



Figure 3. Sketch of the quasilinear diffusion direction and magnitude in vu —v± space 
(i.e. a parallel (z)-perpendicular (x,y) coordinate system). The black contours in the 
x—z plane and the y—z plane are the contours of the gyroaveraged distribution function 
and the brown contours are the characteristics of the operator L(f ). The diffusion 

direction is tangential to the characteristics, v 2 ^ + (y\\ — ^constant. For Landau 
damping, the intensity of the diffusion is determined by the projection of V^/o onto 
the characteristics of the operator L(f ) = at v\\ — w/k\\. The diffusion in velocity 
space results in an average increase of the parallel velocity, Svu . The perpendicular 
momentum transfer Svk ± has the direction of fcj_, but its effect on the distribution 
function vanishes due to the fast gyro-motion (the averaged perpendicular acceleration 
represented by the pink arrows in the yellow circle in the x-y plane vanishes). 




Figure 4. Mechanism for perpendicular momentum transfer by Landau damping of 
the parallel electric field En . The oscillating magnetic field B± is induced by Faraday's 
law. The parallel velocity results in the Lorentz force (F) whose direction is the same 
as kj_, and its magnitude is ZeEu (k±v\\/uj) = ZeE\\(k±/k\\). 



11 



3. Momentum transfer by non-resonant interaction 

The increase in the poloidal mode number is important to determine the location of the 
resonance as shown in Fig.l. In the eikonal limit, the poloidal mode number of the LH 
wave (m ~ kgr in a circular tokamak) is determined by the poloidal variation of the 
determinant of the dispersion relation Dq along a ray path [XT] , 

dm _ dD (uj,m,ne,T e ,B^) / dP (uj, m, n e ,T e , B$) 

dt ~ 06 I dco ' [ } 

Here, in the electrostatic limit, D ~ Sk\ + Pkj, S ~ 1 + u^/Vl 2 . — Upju 2 ~ 0(1) and 
P ~ —Up e /u 2 ^> S are the components of dielectric tensor, and u pe and u P i are electron 
and ion plasma frequency, respectively. Figure 1 shows that as the LH wave propagates 
from the low field side launcher in an inhomogeneous tokamak, the increase in m can be 
as large as ng at the resonance position. Here, n ~ k^R and q ~ are the toroidal 
mode number and the safety factor, respectively. 

As we discussed in Section [2j the poloidal wave momentum is transferred to the electrons 
mainly by resonance, but there is another mechanism that makes the poloidal number 
much larger at the resonance than at the launcher. The origin of the increased poloidal 
momentum is the external force required to keep the density n e , the temperature T e , 
and the static magnetic field constant in time in the dispersion relation. We assume 
that these parameters are fixed in the dispersion relation because the transport and the 
resistive time scale are much longer than the propagation time of the wave. The wave 
exerts a non-resonant force that can affect the evolution of the background profile. In 
general, this non-resonant force is much smaller than the resonant one by a factor of 
Vg/i^fL) < 1 where L is the characteristic length of variation of the background, v g is 
the group velocity of the wave, and 7 is the wave damping rate at the resonance region. 
However, the accumulated momentum transfer by the non-resonant force along the ray 
path is not negligible. 

The nonlinear forces due to the RF wave have been investigated in previous works 
[T2T ITBl EH US HHJ 02] • For example, the nonlinear force exerted by the wave has been 
calculated for a tokamak by neglecting the gradient of the magnetic field compared to 
the gradient of the density and the temperature [15] . In a steady state, these nonlinear 
forces must exactly balance the momentum increase of the wave. Thus, the wave takes 
momentum from the plasma as it propagates due to the plasma inhomogeneity, leading 
to the increase in the wave poloidal momentum. This momentum is given back to the 
plasma at the resonance position. Consequently, the wave has redistributed the poloidal 
momentum of the plasma. The effect that this has on the poloidal rotation is small due 
to the strong poloidal collisional damping in a tokamak [32J. 
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4. Radial particle flux by LH wave 

In this section, we investigate different types of electron radial drifts that can be 
induced by the resonant momentum transfer from the LH wave. The Lorentz force 
and the collisional friction in the parallel direction can balance the toroidal force due 
to waves. The dominant radial electron drift comes from toroidal momentum transfer 
in the perpendicular plane (Sec. 14. ip . The Lorentz force that results from the pinch 
is comparable (0(100%)) to the LH wave momentum source, giving a sizeable radial 
pinch (0(lmm/sec)) that has an outward direction in tokamaks. Other radial pinches 
induced by the wave parallel momentum transfer are relatively small. The passing 
electron pinch caused by the resonance gives a Lorentz force which is 0(10%) of the 
LH wave momentum (Sec. 14.21) . and the Ware-like LH wave induced pinch by trapped 
electrons [25] is associated with only 0(1%) of the LH wave momentum transfer (Sec. 

KB- 



4-1. Outward electron pinch due to perpendicular wave momentum 

The quasilinear term due to the LH wave in the Fokker-Plank equation gives rise to a 
correction to the electron distribution function, F e = f' e — f e , where f' e and f e are the 
electron distribution function with and without LH wave respectively. For convenience, 
we write F e as a function of total energy E = |m e t> 2 — e$, where $ is the background 

2 

potential, magnetic moment p = and the gyro-phase angle a. The equation for 

F e in these variables is 

c)F 3<t> f)F f)P 

~dt~ e ~dt~dE t + ^ V||Fe + Vd ' VFe + ^^cf = Ce{Fe) + Qlh ^> ( 2 °) 
where C e (f e ) is the linearized collision operator to the order of interest (i.e. C e (F e ) = 
C ee (F e ,f e ) + C ee (f e ,F e ) + ]T\ C ei (F e , /j) ) , and v d is the VB and curvature drift. In 
( 120]) we only consider the long wavelength and slowly evolving piece of the distribution 
function because the quasilinear term affects mainly the background distribution 
function. The size of the first and second terms are determined by the gyro-Bohm 
transport time scale, d/dt ~ D g B/a 2 ~ efvu/a ~ y^mi/m e elv te /a, making it much 
smaller than other terms in ( [20]) . Here, D g s = e%PiVi is the gyro-Bohm diffusion 
coefficient, a is the minor radius, is the ion mass, v te and v t % are the electron and ion 
thermal velocities, and e e = p e /a -C 1, = pi/a ~ a/ mi/m e e e 1 are the small ratios 
of electron and ion Larmor radius over the radial scale length, respectively. The third 
term in (120]) is of order v te F e / (qR) where q is the safety factor and R is the major radius. 
The fourth term in ( 1201) is smaller than the third term by (B/Bg)e e . The gyro-motion 
term ^e^f is dominant over other terms (i.e. v\\S7\\F e / (f2 e ^) ~ ae e /{qR) 1 and 
C(F e )/ (fi e ® ~ ^e/fie < !)• Then the lowest order equation is trivial, Vt e ^f = 
(i.e. F e0 = (F e ) a ), and the next order equation is 

U||V||F e0 + Q e — ^ = C e (F e0 ) + QlhUD- (21) 
da 
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Here, we have neglected the time derivative term and the perpendicular drift term. The 
gyro-phase independent part can be obtained by taking the gyro-average of the next 
order equation, 

i>l|V||F e0 = C e (F e0 ) + (Q LH (f e )) a ■ (22) 
The quasilinear term balances with the collision operator and the parallel streaming 
term. The gyro-phase dependent part, F e = F e — (F e ) a , is obtained from the gyro-phase 
dependent contribution to equation (l2~TI) . giving 



n e ^ = Q LH (f e )-(Q LH (f' e )) c 



e 

Its solution is 



-Mr±-<r ± >j + --f + ^-(r l[ -(r ll )j 



v_l <9v_i_ a v± da dv\\ 



(23) 



P el = ^— / da 



AP/ + AP? 



■^(r± - (Tj.>J + --sf + 7 r- (r n -<r„> a ) 



PO ( 24 ) 
n e m e v te ii\ l e 

Thus, the collisional toroidal friction due to the gyro-phase dependent piece of the 
distribution function is much smaller than the corresponding RF force, 

f dv\m e Rv^)C{Ki) ~ (AP/ + AP;) , (25) 

and most of the perpendicular momentum transfer is balanced by the Lorentz force 
/ dv^m.Rv^n^ from Q23&. 

The radial particle flux can be obtained from 

(T e • W) s ~ l^J dv 3 F e0 v d • + l^J dv 3 F\ lV± p ■ , (26) 

where ip is the poloidal magnetic flux, (...) s is the flux-surface average (see |Appendix ) 



and T e is the electron particle flux due to the correction F e . From the steady state 
Fokker-Planck equation given in (j2ip . taking the moment (m e f</,P) and a flux-surface 
average of the resulting moment equation, we can relate the radial pinch (r e • VV>) S to 
the correction F e by 

c \ c 

+ ({i ■ f>)R J d\m e v\\ [C(P e0 ) + (Q LH (f' e ))Sj . (27) 




To obtain this equation (127|). we use that the first and second term on the left hand side 
of (|2ip give the first and second term on the right hand side of (l26j) . respectively. The 
right hand side of ( 1271) is obtained by decomposing the right hand side of ( 12~TT) into the 
gyro-phase dependent and gyro-phase independent pieces, 

C e (P e ) + Q LH (fe) * [Qim{& ~ (Qlh(K)U + + (QlhW))*] ' ( 28 ) 
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The second term on the right hand side of (127)) is the parallel force balance from the sec- 
ond term on the right hand side of (1281) . which will be discussed in the next subsection. 
The first term on the right hand side of ( )27j) is the toroidal projection of the perpendic- 
ular wave momentum transfer, (AP^- + AP?) , which comes from the first term on the 
right hand side of (128)) . In ( 1281) we have already neglected the perpendicular collisional 
friction (see ( 125)) ). 

The collisions transfer most of the parallel wave momentum to the ions, but the rest of 
the toroidal angular momentum (e.g. n^R — n\\(B ( / ) /B)R ~ 0.9 in Fig. 1) remains and 
it has the opposite toroidal direction to the original toroidal angular wave momentum, 
giving an electron outward pinch that is opposite to the inward pinch predicted in pre- 
vious works [2} [25], [26] . Physically, the outward radial pinch comes from the effect of 
the perpendicular wave momentum transfer AP^ + AP^ on the gyro-motion (see Fig. 
5- (a)). This electron pinch is still very small compared to the Ware pinch [27]. For 
example, if 1MW of LH wave power is locally absorbed in a volume of 0.1m 3 where the 
plasma density is 10 20 m -3 , the poloidal magnetic field is Bq = 0.5T, and the refractive 
index is — n\\(B^/B) = 1, then the electron outward radial pinch is about 4 mm/s 
which is a hundred times smaller than the Ware pinch. 



Jf.,2. Passing electron pinch due to parallel wave momentum 

To solve for the gyro-phase independent perturbation (F e ) a due to the LH wave in ( 1221) . 
we use a subsidiary expansion of F e0 = P e ° + F^ + in the small ratio of the collision 
frequency over the transit frequency in the banana regime. The lowest order equation 
is V||P e ° = 0, implying that P e ° is a flux function. The next order equation is 

v^\\Fl, = C e {F! ) + (Q LH {f e )) a . (29) 



Taking a bounce average of ( 1291) (see Appendix C), the left hand side of ( E291) vanishes 



(C e (F° ) + Q LH (f>)) b = 0. (30) 
According to ( jC.lj) . this equation is equivalent to 

f (Ce(F! )+Q LH (f e ))\ =0. (31) 

In general, the solution P e ° to equation ()3"Tj) does not make the second term on the right 
hand side of ( 127)) vanish, giving a non-zero radial pinch due to the passing electrons. 
This imbalance comes from the variation of v», B and R along the orbit, which is of the 
order of the local aspect ratio, 0(r/R). For the electrons resonant with the LH wave, 
the effect of the change of vn along the orbit is negligible because most resonant elec- 
trons have much larger parallel velocity than perpendicular velocity (i.e. small magnetic 
moment, fi ~ 0) due to the high phase velocity of the wave in the parallel direction. 
The non-vanishing contribution to the radial pinch is due the competition between the 
localized wave power absorption within a flux surface and the collisions that occur over 
the whole flux surface. 



15 



Physically, this pinch can be explained by how the passing orbit of a single elec- 
tron is changed by the resonance. The canonical angular momentum of the electron 
ip* = if) + Iv\\/Q, e determines the radial deviation of the electron orbit from the flux sur- 
face ip due to the curvature and V-B drifts. Here I = RB^ is a flux function to lowest 
order. After the resonance with the negative k\\ of the LH wave, the absolute value of the 
negative velocity of the resonant electrons is increased by \Av\\ | due to the absorbed wave 
power. Accordingly, the change of the canonical momentum is Aip* = IAv»/fl e < 0, 
where the gyrofrequency Q e is evaluated at the local resonance point within the flux 
surface. Assuming the low frequency collisions cause the resonant electron to lose its 
momentum only after many transits, we can use the temporally averaged radial location 
to describe its radial motion. The increase in the transit averaged (bounce averaged) 
radial position of the particle is 

AW ' = <* > '-<^ = A *'-{(®),) t -<@) 1 >} (32) 

IAv\\ m e cI(§ 2 dl/B § x dl/B , 



Vt e e \ r b2 r bl 

- ,Av * (34) 

Here, the values with the subscripts 1 and 2 are before and after the resonance, respec- 
tively, Tf, is the bounce time, and from ( 152"]) to (1331) the equation (IC.ll) is used. From (1331) 
to ( |34l) . we neglect the radial displacement due to the poloidal variation of the parallel 
velocity because of the small magnetic moment (/x ~ 0). The flux surface averaged value 
in (1331) is approximated by that at the magnetic axis using a small inverse aspect ratio 
expansion. The frequency Q e o is the gyrofrequency at the magnetic axis. Equation 
( |34l) means that the temporally averaged particle radial flux due to the resonance is 
negative for a low field side resonance (inward radial pinch) and positive for a high field 
side resonance (outward radial pinch). It is shown in Fig. 5-(b). The increase in the 
curvature drift due to the increase in the parallel velocity after the resonance results 
in the different passing orbits depending on the resonance location on the flux surface. 
This radial drift is included in the second term on the right hand side of (1271) as the 
competition between the localized wave power absorption within a flux surface and the 
collisions that occur over the whole flux surface. For a typical small inverse aspect ratio 
tokamak, this imbalance is small, about 10% (0(r/R)) of the total momentum transfer. 



4-3. Trapped electron pinch due to parallel wave momentum 

For trapped electrons, since odd functions in vn vanish under the bounce averaging 
according to (1C.2I) . equation (1301) becomes 

(C e (F° )) b = -(QiT(f' e )) b . (35) 
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The distribution function F e0 is an even function of v», because the bounce averaged 
quasilinear term is even. Two trapped electrons at the outer-midplane that have oppo- 
site parallel velocities have the same electron Landau damping resonance at the same 
local point in their banana orbits due to the small electron banana width (see Fig. 5- 
(c)). The non-zero F e0 due to the LH wave can be understood as follows: a trapped 
electron is accelerated only when it resonates with the wave, that is, when its velocity 
is the same as the wave phase velocity, and this acceleration continues every transit 
until it collides. As a result, it has an open trajectory that moves inward every bounce. 
There is no net gain of toroidal angular momentum for the trapped electron because F e0 
is an even function, but there is a gain of canonical angular momentum that leads to a 
Ware-like LH induced pinch. This pinch is the contribution of the trapped electrons to 
the second term on the right hand side of ( l27|) . The LH wave trapped electron pinch 
is tiny, because the power absorption by trapped electron is less than 1% of total wave 
power due to the small size of the population of trapped electrons resonant with the LH 
wave phase velocity (co/ku ~ 3v t h — I0v t h) . It results in a very small contribution to 
the radial pinch (approximately less than O.lmm/s). 

The mechanism behind the radial pinch by trapped electrons is similar to the mech- 
anism of the outward radial pinch due to the gyro-motion described in section 14.11 
Instead of considering the effect of the acceleration on the gyro-motion, one needs to 
consider its effect on the banana orbit (compare Fig. 5- (a) and Fig 5-(c)). 

5. Conclusion 

In this paper, we have proven that wave-particle momentum transfer by resonance 
happens through both the parallel motion and the gyro-motion. Only considering the 
parallel motion leads to incorrect results when evaluating toroidal angular momentum 
transfer. The toroidal momentum carried by the parallel motion is rapidly dissipated 
by collision with ions. The perpendicular force is balanced by an electron radial pinch 
rather than collisions. For the LH wave in tokamaks, the difference between the toroidal 
wave vector and the toroidal component of the parallel wave vector, represented by 
n^R — n\\(B^/B)R, determines the radial pinch of the electrons. Typically, the high 
poloidal wave number at the electron Landau resonance induces a pinch with the 
opposite sign to the pinch that one would have expected for a wave with no poloidal wave 
number and the same toroidal wave number. For counter-current direction LH wave 
momentum input, this pinch is 0(lmm/sec) and it generates an additional outward 
radial electric field that makes the flux of electrons and ions ambipolar. This radial 
electric field gives an E x B ion flow with the opposite direction to the momentum 
source. Eventually, after receiving the parallel momentum by collisions with electrons, 
the ion velocity will acquire the direction of the momentum source. The ions achieve the 
toroidal angular momentum that the LH wave contained originally through two main 
channels: ion-electron collision for the parallel direction motion, and the Lorentz force 
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Figure 5. (a) Sketch of the outward radial drift during the gyro-motion due to 
the perpendicular momentum transfer, (AP^ + AP^ . The direction of the toroidal 
component of the perpendicular force is the same as the direction of the plasma current 
Jtjj. (b) Sketch of the poloidal cross section of the trajectories for two passing electrons 
that receive wave parallel momentum at the inner-midplane (blue) and the outer- 
midplane (red), respectively. The dashed black line is the transit orbit before the 
resonances. Due to the increased curvature drift (upward direction) after the resonance, 
the passing electron orbits are different depending on the resonance location. The 
temporally averaged radial flux is outward for a resonance at R < Rq and inward for 
a resonance at R > i?o- (c) Sketch of the poloidal cross section of the trajectories for 
trapped electrons that receive wave parallel momentum. The dashed black line is a 
flux surface. For both signs of v\\ at the outer-midplane on the flux surface, the orbits 
move inward and increase its width with every bounce. 



due to an outward radial ion pinch following the electron pinch. 



Appendix A. Wave power absorption by quasilinear diffusion 

To evaluate the wave power absorption in fl9]), we utilize the Bessel function expansion 
for the sinusoid phase, 

= 5>^J n (A) (A.l) 



sm^= -J>^j;(A) (A.2) 

n 

cosrje lXsinri = V^e in V n (A), (A.3) 



and the sifting property of the phase average 

»2tt p2n 



/ dr]e -^(v+nr)-sin( v )) = / ^V" e -^+^)j / ( A )^ e ^j n (A) (A.4) 
Jo Jo I 

= 2nJ2e- mQT J 2 n (X). (A.5) 
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Then, the power absorption P a bs,\a for a single mode can be evaluated [5J EJJ as 

x poo r°o / mv 2 

Pabs = Pabsfr = J dv\\ J dv±27rv_L(—-Q(f) 

(A.6) 

giving 

'2„2 poo 



Pabs,k = — dv \\ dv L 27cv L ^——L(vl5(u-k l \v l \-ntt)\xk,n\ 2 L{fo))(A.7) 

m J-oo J° n 

2„2 poo poo 

/ dv\\ dv±2nv ± ^mv 2 ± 5(uj - k\\v\\ - nVL)\xk,n\ 2 L(fo)- (A. 8) 



7rZ e 



Appendix B. Wave momentum transfer by resonances 

The toroidal momentum transfer in the parallel direction ( ITT]) can be written by inserting 
the flux T|| §6§ and the perturbed fluctuated distribution function / k (jSJ). 

/oo poo p2n 

dv\\ I dv±v± / da(z ■ 0) 
K oo Jo Jo 

x{B;, || (l-^icos(a-/ 3 )) + B t± fe}/ k 

Z2 2 p poo poo poo p2it 



m 

k 



/OO POO POO PZ7V 

dv\\ I dv ± 2irv ± dTe i{ "- k w v w )T I d<n{z-$) 
-oo Jo Jo Jo 



X 

n 



x ^ e"*^) (jSicii^JJ + y((^k I+ + ^k,-)f/ - E\\V) + (E k>+ - E^)UJ^ . 

The phase i(u — fc||t>|| — nf2)r is averaged out in the r integration except where 
u> — k\\V\\ — nQ = 0. Using the Dirac-delta function to express this resonance condition, 
the momentum transfer becomes 

p\ = ~ ~ — ~~~ / dv\\ f dv±2irv^ 8(cu — — nft)(z ■ 0) 

x ((l - ^) «,„ + =^(* k , + + + J^(E i+ - E U ) 

x [E ki ^-J n + ^((B k ,+ + £*,_)£/ - B k .,V) + (B k , + - B k ,-)(7^J (B.l) 

/oo /*oo 
/ df_|_27rfj_ 
■oo Jo 



irZ 2 e 2 R 



m 

k 



j 2 

x^(w- Vll " • $) J ^\XKn\ 2 Hfo) (B.2) 

^— ' a; 

n 

—Pab.,vR@ -4>)=y\ —PabsMRi* ■ 4>) , (B.3) 

^— ' a; c 

k k 
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where the resonance condition uj — kuVu —nVt = 0, the Bessel function identities nJ n /X = 
(J„+i + J„_i)/2 and J' n = (J n _ a - J n+1 )/2, and Xk,n = ^k,||^n^ + #k,+</n-i + E k -J n +i 
are used from (jRH to (IB~2l . 



The rest of the toroidal momentum transfer in (JT6J) can be obtained by inserting the 
perturbed fluctuated distribution function /k ©• Before doing so, we rewrite f|T6|) as 

/oo /»oo ^"2^ 7 

/ ^ ± 2tt^ / — 
K -00 ^0 ./o l7[ 

k±v\\ 



El 



uj 

k»v\\ 

UJ 



(cos /3(x ■ 0) + sin/3(y ■ </>)) 

+ E* K _) cos /3 - i{El + - E* K _) sin /?) (x • 0) 



+ (K , + K _) sin /3 + i(E* - E*) cos/3 (y • <f>) } + i 



\ UJ 



(Et 



El 



x j (sin 77 cos /3 + cosrisin/3)(x • cj)) — (cos 77 cos /3 — sinr7sin/3)(y • 0)| /k- (B.4) 
Using ([H]) for /k and the Dirac-delta function for the resonance condition gives 

ap; + ap; 

irZ 2 e 2 R 



rn 



E 



x "\ JnE k i 



1 - 





fcllf II 



dv±2irv±5(uj — kuvu — (cos /3(x • 0) + sin /3(y • 0)) 



x 



UJ \ UJ 

i(sm /9(x • 4>) — cos /3(y • cj>))) 
fciiuii n£> 



7 l Z?* I Z7 1 * \ 1 7' ^-L^-L / 771* 771* 



1 - 



UJ 



nZ 2 e 2 R 



UJ 

"OO 



El 



XKnV±L{f ) 



(B.5) 



/oo /*oo 
(ifn / (if j_27rf j_ 5(uj — k\\v\\ — n^) 
-00 Jo 



x (^cos /3(x • cj)) + sin /3(y • cj)) 
k .0 



|Xk,n| 2 ^(/o) 



(B.6) 
(B.7) 

k k 

From step (IB.5P to flB.6p . the resonance condition and the Bessel function identities 
nJ n /\ = (J n+ i + J n _i)/2 and = (J n _i - J„+i)/2 are used. 



y~] — Pabs,kR = Y] —Pabs,k R (cos /3 (it ■ 0) + sin/3(y ■ 0) 



Appendix C. Bounce averaging and flux surface averaging 

The bounce- averag is defined as (X)b — A § ^Ol = i. ^ A|A^ where 6* is the poloidal 
angle and r b (£,fi) = § — ^ ^ is the bounce time. The flux surface average is 

(A), = ±;§ §m = tJ^i because the Jacobian is J = V( ^. w = and the 
normalization factor is r s = j> ^. There is a relation between the flux surface averaged 
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(C.l) 



For trapped particles, assuming that the flux surface average is taken between the 
turning points of the orbit, the bounce averaged value is defined as 



where r s i and r h \ stand for the integration and the bounce time from one turning point 
to the other turning point, and the summation over a = v»/\vu | indicates that the values 
of X for the two parallel velocity signs must be added. Bounce averaging annihilates 
the operator v yVyX for both passing and trapped particles. 
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